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INTRODUCTION 


The  subject  of  thermal  radiation  from,  or  emisaivity  of,  micrometer-sized 
particulate  matter  is  of  general  interest,  but  the  studies  of  this  radiation 
have  met  with  difficulty  in  obtaining  experimental  data  with  unambiguous 
interpretation.  Most  situations  in  which  particles  are  generated  at  high 
temperatures,  such  as  rocket  or  engine  exhausts,  are  quite  complex;  the 
observed  radiation  spectra  are  not  sufficiently  structured  to  specify  the  hoot 
of  variables  which  enter  an  a  priori  calculation  of  particle  emission. 
Typical  emission  parameters  which  may  be  difficult  to  measure  in  combustive 
flows  producing  or  including  particles  would  be:  (a)  particle  size 
distribution,  (b)  particle  composition  distribution,  (c)  particle  structure 
and  shape  distribution,  (d)  corresponding  bulk  indices  of  refraction, 
(e)  spatial  and  temporal  changes  in  these  distributions  and  al3o  in 
temperature,  and  (f)  other  molecular  (gas)  species  which  may  emit  continuous 
or  quasi-continuous  bands  of  infrared  radiation. 

Because  there  are  so  many  variables,  the  general  radiation  charac¬ 
teristics  of  particulate  matter  have  not  been  verified  experimentally.  In 
addition,  there  are  fundamental  theoretical  questions  still  unanswered.  What 

theoretical  model  should  be  used  in  predicting  particle  radiation?*  What  is 
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the  size  limitation  of  a  valid  "classical"  description?  What  is  the  effect 
of  exceeding  that  size  limitation?  This  report  presents  results  of  an  experi¬ 
ment  in  which  the  variables  are  either  measured  or  minimized  and  compares  the 
observed  spectra  with  calculations  based  on  Mie  Theory,  Planck  blackbody 
radiation,  and  bulk  values  of  the  index  of  refraction. 


1.  McGregor,  W.  K.,  "On  the  Radiation  from  Small  Particles,"  J.  Quant. 
Spectropscopy  Radiation  Transfer  19,  659  (1978). 

2.  Kawabata,  A.  and  Kubo,  R.  0.,  Phys.  Soc.  Japan  21,  1765  (1966). 


SECTION  1 
OBJECTIVES 


The  principal  objective  of  this  program  was  to  explore  possibilities  for 
more  refined,  nonobtrusive  diagnostic  techniques  spplicable  to  combustive 
f lows  such  as  a  solid  propellant  rocket  motor  plume.  One  of  the  areas 
investigated  was  the  thermal  emission  from  two-phase  flows.  The  experiment 
was  designed  to  exert  as  much  control  as  possible  over  the  various  parameters 
which  affect  the  thermal  radiation  spectra  emitted  by  suspended  particulate 
material.  Therefore,  the  primary  objective  was  to  measure  accurately  the 
thermal  radiation  from  a  well-characterized  two-phase  flow  and  then,  if 
possible,  to  use  these  data  to  refine  theoretical  formulations  for  predicting 
infrared  emission  from  such  flows. 
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SECTION  2 

EXPERIMENTAL  APPROACH 


2.1  TWO-PHASE  FLOW  APPARATUS 

The  experimental  approach  i8  in  principle  8imple  and  straightforward. 
The  greatest  difficulty  is  in  constructing  a  gaa  and  particle,  or  two-phase, 
flow  apparatus  with  a  constant  flow  velocity  in  which  the  composition,  size, 
temperature,  and  particle  density  are  controlled.  Figure  1  is  a  schematic 
diagram  of  the  two-phase  flow  apparatus  which  has  evolved.  The  basic  com¬ 
ponents  of  the  system  are  a  powder  hopper  with  an  auger  feed,  a  particle 
suspension  tank,  aid  a  gas  heater.  The  particle  feeder  is  a  commercial  unit 
which  can  hold  30  4  of  powdered  material.  The  feeder  transports  the  particles 
with  a  variable  speed  auger  screw.  A  gas  heater  is  constructed  out  of  roughly 
A  m  of  12.5-mm  diameter  stainless  steel  tubing  bent  into  a  coil.  Short  pieces 
of  smaller  diameter  stainless  steel  tubing  were  inserted  before  the  coil  was 
formed  to  improve  heat  transfer  to  the  gaa  by  breaking  up  laminar  flow  and 
providing  greater  surface  area.  Cables  to  an  arc  welder  capable  of  providing 
up  to  200  A  are  connected  to  either  side  of  the  coil  which  is  isolated  from 
electrical  ground  by  ceramic  feed-throughs.  The  coil  is  then  enclosed  in  an 
aluminum  box  which  is  thermally  insulated  with  zirconium  oxide.  When  the  tube 
is  electrically  heated,  the  N£  gas  flowing  through  the  tubing  can  reach 
temperatures  in  excess  of  1,300  K.  The  main  gas  stream  i3  introduced  into  the 
system  at  about  180  kPa  (25  psi).  A  smaller  N^  gas  flow  is  maintained  through 
the  holding-suspension  tank  to  transport  the  particles  through  the  exit 
nozzle.  The  nozzle  throat  was  10  mm  in  diameter,  and  pressure  behind  it  was 
principally  due  to  the  carrier  gas  flow.  This  overpressure  was  typically  held 
between  10  and  20  kPa,  since  the  particle  feeder  was  not  originally  designed 
to  hold  pressure.  Calculations  based  on  this  pressure  and  measurements  made 
with  anemometers  indicate  normal  flow  velocities  of  50  m/s  at  the  exit  plane. 


2.2  OPTICAL  SYSTEM 

Once  a  flow  of  heated  particles  is  achieved,  it  is  very  straightforward 
to  set  up  instrumentation  for  the  detection  and  spectral  analysis  of  infrared 
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Figure  1.  Two-Phase  (Gas/Particle)  Flow  Generator. 


radiation.  Since  nitrogen  gas  does  not  have  infrared-active  modes,  any 
radiation  observed  must  be  due  to  the  particles.  The  optical  path  consists  of 
a  collecting  lens  which  focuses  the  emitted  radiation  onto  the  entrance  slit 
of  a  1/4-m  scanning  monochromator,  and  a  curved  mirror  which  focuses  the 
radiation  from  the  exit  slit  onto  an  InSb  detector.  A  chopping  wheel  is 
placed  in  front  of  the  monochromator.  Output  of  the  detector  and  preamplifier 
is  measured  with  a  lock-in  amplifier  synchronized  to  the  chopping  frequency. 

Figure  2  shows  photographs  of  the  apparatus.  In  Figure  2(a),  the  par¬ 
ticle  feeder  is  located  on  the  left.  On  the  right  is  an  instrument  rack  con¬ 
taining  the  monochromator  and  associated  electronics.  In  the  lower  middle 
section  of  the  picture  is  the  ga3  heating  coil  with  the  insulation  box 
removed.  Behind  the  coil  is  the  suspension-holding  tank  with  the  exit  nozzle 
on  top.  Figure  2(b)  is  a  different  perspective  of  the  same  apparatus,  showing 
the  holding-suspension  tank  in  the  lower  middle.  The  heating  coil  is  on  the 
right,  with  the  insulation  box  now  in  place. 

In  addition  to  the  infrared  radiation  spectral  measurement,  other  optical 
systems  monitor  relative  changes  in  particle  density  in  the  flow  stream. 
Scattering  of  HeNe  laser  light  by  the  particles  is  detected  by  a  photo¬ 
multiplier  tube  and  lens,  and  the  transmitted  beam  is  monitored  by  a 
photodiode.  Figure  3  is  a  schematic  diagram  of  the  optical  arrangement  for 
making  all  the  various  measurements,  although  they  were  generally  not  carried 
out  simultaneously.  All  the  optical  diagnostic  measurements  are  centered  at  a 
height  of  10  mm  (1  nozzle  diameter)  from  the  nozzle  exit  plane.  With  the 
laser  unfocused,  photon  correlation  spectroscopy  techniques  could  be  used,  and 
with  a  focused  beam,  discrete  particles  could  be  counted.  A  blackbody 

calibration  source  i3  used  to  determine  the  absolute  instrument  response  of 
the  monochromator  and  detector  (Fig.  3).  Due  to  the  effects  of  second  order 

diffraction,  it  is  necessary  to  use  cutoff  filters  and  measure  the  spectrum  in 

three  wavelength  segments:  1.5  to  2.4  ym,  2.4  to  3.6  urn,  and  3.6  to  5.5  ym. 
The  long  wavelength  cutoff  is  due  to  the  cooled  InSb  detector  response. 

The  laser  light  scattering  signal  from  the  photomultiplier  tube  is 
amplified  and  recorded  on  FM  tape  with  a  frequency  response  of  0-10  kHz.  For 
the  photon  correlation  spectroscopy  measurements  the  signal  is  played  back 

into  an  analog-to-digital  fast  Fourier  transform  computer  (SD-360),  and  the 
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Figtir;  2(a  and  h). 


^holographs  From  Different  Perspectives  of  the  Two- 
Phase  Tlow  Generator  and  the  Infrared  Radiation  Detect  in 
Instrumentation.  (a)  Detection  Instrumentation  is  right 
with  Flow  Generator  left.  (b)  Detection  Instrumentation 
is  behind  Flow  Apparatus  in  foreground. 
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Schematic  Diagram  of  Optical  Measurement  Geometry  Looking  Along  the  Axis  of  the 
Two-Phase  Flow. 


results  output  to  an  x-y  plotter.  For  the  discrete  particle  detection  the  FM 
tape  is  speed-scaled  down  and  played  back  into  an  optical  chart  recorder.  The 
laser  attenuation  is  monitored  by  a  photodiode/op-amp  detector  input  directly 
to  a  strip  chart  recorder. 

2.3  PARTICULATE  MATERIAL 

Only  one  type  of  particulate  material  has  been  used  in  the  flow  apparatus 
-  a  pure  carbon  powder  available  commercially.^  The  size  distribution  of  this 
material  aa  single  particles  was  measured  on  a  scanning  electron  microscope, 
as  shown  in  Figure  4.  Nearly  all  the  powder  has  an  initial  size  of  less  than 
0.5  urn  in  diameter.  There  was  the  question  of  whether  this  remains  an 
accurate  description  of  particle  sizes  in  the  gas  flow  since  agglomeration 
could  take  place  during  transport  and  mixing  with  the  nitrogen.  Samples  were 
taken  from  the  two-phase  flow,  by  passing  a  clean  tantalum  surface  through  the 
flow  stream.  These  tantalum  discs  were  examined  with  a  acanning  electron 
microscope,  and  while  individual  carbon  spheres  on  the  order  of  0.5  urn  were  in 
evidence,  there  were  also  a  significant  number  of  large  agglomerated  particles 
up  to  15  urn  in  diameter.  Figure  5  (a  and  b)  3hows  representative 
photomicrographs  taken  with  the  scanning  electron  microscope. 

Some  indication  of  the  degree  of  this  agglomeration  was  needed,  since  the 
spectral  data  would  be  a  function  of  the  particle  size  distribution.  The 
simplest  technique  for  determining  average  size  i3  to  measure  simultaneously 
laser  beam  attenuation  and  scattering.  By  focusing  the  laser  beam  to  a  small 
enough  diameter,  it  is  possible  to  use  a  photomultiplier  tube  to  detect  scat¬ 
tered  laser  light  from  a  single  particle  rather  than  a  group  of  particles.  By 
counting  the  average  number  of  discrete  particles  under  normal  flow  conditions 

and  calculating  the  scattering  volume,  the  density  of  particles  was  directly 

4  3 

measured  to  be  ~5  x  10  /cm  .  This  number  can  be  coupled  with  the  attenuation 


3.  Thermax  Powder  from  R.  T.  Vanderbilt  Co.  Inc.,  6273  E.  Slausen, 
Los  Angeles,  CA. 
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Figure  4.  Distribution  of  Single  Carbon  Particle  Diameters  Used  in  the  Two-Phase  Flow. 


Figure  5(a  and  b), 


Photographs  From  electron  Microscope  of  Particulate 
Samples  Collected  in  the  Two-Phase  Flow  at  low  Flow 
Rates.  In  5a  the  scale  bar  is  5  um.  Tn  5b  the  scale 
bar  is  1  ym  and  the  individual  particles  can  be  seen 
as  0.5  urn  or  less  in  diameter.  The  question  is 
whether  such  agglomerates  are  present  in  the  flow. 
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of  the  laser  beam  as  it  passes  through  the  particle  flow  by  the  equation: 


-ln(I/I  )  =  p  •  o  •  L  (1) 
o 

where  I  and  I  are  the  initial  and  attenuated  laser  intensities,  L  is  the 
path length  through  the  two-phase  flow,  and  p  and  a  are  the  average  density  and 
extinction  cross  section  for  the  particles.  Average  attenuation  of  the  laser 
is  2.5  +  0.5  percent.  Using  5  x  10  /cm  for  p,  a  can  be  calculated  to  provide 
an  average  particle  diameter  of  2.6  pm. 

A  more  sophisticated  technique  was  also  employed  to  obtain  an  average 
particle  size.  By  recording  the  Rayleigh  (elastic)  scattering  from  a  group  or 
population  of  particles,  correlation  spectroscopy  can  be  used  to  determine 
their  diffusion  rate.  Detailed  explanations  of  this  technique  can  be  found 
elsewhere. Only  a  concise  outline  is  given  here.  Consider  the  scattered 
electric  field  from  a  single  particle.  From  Figure  6,  £0  is  the  initial  light 
wave  vector  and  kg  the  wave  scattered  from  the  particle.  The  electric 
field  due  to  scattering  from  this  particle  can  be  written  as 


£■  =  £  exp  i[£  •  r  .+  i<  •  (ft  -  r .)  -  wt]  (2) 

j  orLojs  o  j 

where  fj  is  the  position  vector  from  the  origin  to  the  jth  particle,  and  is 
the  position  vector  to  the  detector.  A  simplification  is  made  by  defining  a 
scattering  wave  vector,  ft,  by 

|*|=  SQ-  S'  (3) 


4.  Renner,  S.  S.,  Bernard,  0.  M. ,  and  Jerskey,  T.,  Acta  Astronautica  2»  93 
(1976). 

5.  Clark,  N.  A.,  Lunacek,  J.  H.,  and  Bendek,  G.  B.,  Am.  J.  Phys  38^  575 
(1970). 

Schaefer,  D.  W.  and  Berne,  B.  0.,  Phys.  Rev.  Letters  28,  475  (1972). 
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From  Figure  6  it  can  be  shown  that 


i 


|i?|  =  2  k  sin(6/2) 


Substituting  Eq.  (3)  into  Eq.  (2)  yields 


fj*  e*P  U  $  *  fj)} 


t 

where  E3  =  EQ  exp  {i(kg  •  R0-  ut)}  .  The  total  electric  field  at  the  detector 
is  merely  the  sum  of  the  scattering  fields  from  N  individual  particles,  or, 


V  f  j  =  K  £  e*e  [i(*  *  vj 


A  useful  expression  would  link  the  electric  field  to  the  power  spectrum  of  the 
intensity  of  light  at  the  detector.  A  connection  is  found  in  the 
Wiener-Khinchine  theorem  which  states  that  the  power  spectrum  is  the  Fourier 
transform  of  the  autocorrelation  function. 

The  autocorrelation  function  measures  the  degree  of  coherence  a  function 
has  with  a  copy  of  itself  shifted  in  time,  by  an  amount  t.  If  "n"  is  a  par¬ 
ticle  at  time  t,  and  "m"  a  particle  at  time  t  +  t,  the  field  correlation  func¬ 
tion  is  given  by 


N  N 

/etU)  Et*  (t«)dt  =  2  S/EsEsexP  j*[*  '  "n(t)-rm(UT)]Jdt  (7) 

n=l  m=l  ’ 

Because  the  time  dependence  is  evident  only  in  the  change  in  particle 
position,  the  integral  can  be  changed  from  a  time -dependent  to  a  position- 
dependent  integral,  and  Eq.  (7)  rewritten  as 


v  ■  r 


(8) 


N  N 

£  SA  E8exp  I1*  •pn(t)“r™(t+T)]  [  dt 

n=l  m=I 


=  NE2  ^  ^  exp(i£  •  Ar^ 

J 

where  Arj  =  if  j ( t )  -  rj(t+x).  The  ensemble  average  over  Arj  is  given  by 

^  j  exp  (i$  *  A^)  =  J"  f(Ar,x)  exp  ^i(£  *  Arj)J  dAr  (9) 

where  f(Ar,x)  is  the  conditional  probability  that  a  particle  moves  A?  in  a 
time  period  t.  In  the  model,  the  only  change  in  position  through  time  t  is 
due  to  random  Brownian  motion;  thus  the  probability  of  a  particle  moving  Ar  in 
a  time  x  is  given  by^ 


f ( Ar, t)  =  (4*Dx)'1/2  exp(~  |at  |2/AOt)  (10) 

Where  D  is  the  diffusion  coefficient  for  Stoke 'a  flow  and  where  molecular  mean 
free  path  is  small  compared  to  particle  size, 


D  =  kT/3irud 


(ID 


where 

k  is  the  Boltzmann  constant 
T  is  the  temperature 
u  is  the  coefficient  of  viscosity 
d  is  the  average  particle  diameter. 


7.  Present,  R.  D.  Kinetic  Theory  of  Gases  (McGraw-Hill,  New  York,  1958), 
p.  158  . 
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Combining  Lq's.  (7)  through  (11)  and  integrating  gives  the  field  correlation 
function  as 


jl E(t)E*(t+r)dt  =  NE2  exp(-K2Dr) 


(12) 


Using  the  Weiner-Khinchine  theorem,  the  Fourier  transform  of  Eq.  (12)  will 
provide  an  expression  for  the  power  spectrum, 


P  (u) 
s 


=  2 71  / 


NE2  exp(-K2Dx)  exp(-iwx)dx  = 


NE2K2D/ir 
3 _ 


(K2D)2 


(u>)' 


(13) 


This  expression  for  the  power  spectrum  is  conveniently  in  the  form  of  a 
Lorentzian  curve  given  by 


Y  =  - 5 - T  (14) 

(HWr  +  * 

where  HW  is  the  "ha  If  width"  of  the  curve.  So  far  the  derivation  has  dealt 
with  the  field  correlation  function  for  simplicity.  However,  the  measurable 
quantity  is  the  scattered  light  intensity,  not  the  field  strength.  Hence,  the 
field  correlation  function  (Eq.  12)  is  squared  to  obtain  the  intensity  corre¬ 
lation  function.  Since  squaring  an  exponential  function  simply  multiplies  the 
argument  by  a  factor  of  two,  the  form  of  the  power  spectrum  is  preserved 
(neglecting  dc  terms).  In  going  from  the  field  to  the  intensity,  simply 
substitute  2K^d  for  K^D  in  Eq.  (13).  This  implies  that  the  intensity  power 
spectrum  is  twice  as  wide  as  the  field  power  spectrum.  This  derivation  makes 
the  necessary  theoretical  connection  between  the  intensity  power  spectrum  and 
physical  parameters.  Operationally,  the  light  intensity  power  spectrum  is 
obtained  by  taking  the  absolute  value  squared  of  the  Fourier  transform  of  the 
measured  scattered  light  intensity  as  a  function  of  time.  Therefore,  the 
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power  spectrum  of  the  scattered  light  intensity  as  a  function  of  time  yields  a 
Lorentzian  curve  with  a  halfwidth  given  by: 


7  4k^  sin^(6/2)kT 

HW  =  2ICD  =  — - - = -  rad/sec  (15) 

pdir 


figure  7  is  a  graph  of  the  scattered  light  signal  (time  series)  from  the 
photomultiplier,  and  the  corresponding  power  spectrum  from  which  the  halfwidth 
of  50  Hz  yields  an  average  diameter  of  3.4  pm  through  Eq.  (15).  This  tech¬ 
nique  was  perfected  and  calibrated  by  first  performing  a  similar  measurement 
on  polystyrene  latex  spheres  suspended  in  water.  Since  the  latex  spheres  were 
commercially  prepared  and  measured  to  be  a  certain  diameter  (0.312pm  in  this 
case),  the  degree  of  absolute  error  in  the  measurement  could  be  assessed  and 
was  found  to  be  ~7  percent.  However,  all  the  usable  data  obtained  on  the  two- 
phase  flow  apparatus  were  taken  at  low  flow  rates.  Velocities  on  the  order  of 
2  cm/s  were  indicated  from  the  widths  of  single  particle-scattering  pulses. 
When  normal  flow  conditions  were  established  (~50  m/s),  a  large  spike  centered 
at  0  Hz  in  the  power  spectrum  buried  the  data.  Therefore,  utilization  of 
these  data  requires  the  assumption  that  the  size  distribution  is  independent 
of  the  flow  conditions.  The  indicated  particle  3ize  of  3.4  pm  is  therefore 
considered  to  be  in  reasonably  good  agreement  with  the  previous  measurement  of 
2.6  pm,  taken  under  nprmal  flow  conditions  independently. 

Some  assessment  of  the  degree  of  agglomeration  is  achieved  on  the  basis 
of  these  in  situ  measurements  and  calculations.  However,  there  is  still  a 
certain  degree  of  arbitrariness  in  specifying  the  particle  size  distribution, 
since  the  measurements  reflect  only  on  some  sort  of  "mean"  size.  In  the 
absence  of  more  information  about  the  particle  size  distribution  in  the  flow, 
the  analysis  of  the  data  presumed  a  standard  two-parameter,  log-normal  distri¬ 
bution.  One  parameter  is  obtained  from  the  particle  size  measurements,  which 
are  equated  with  the  geometric  mean  of  the  distribution  since  they  are  based 
on  light  scattering  pnd  extinction,  and  therefore  weigh  the  results  by  the 
particle  cross-section  (proportional  to  d^  in  this  case  where  d  ~  X).  The 
second  parameter  is  derived  from  the  hypothesis  that  the  most  numerous 
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Figure  7 


Scattered  Laser  Light  Signal  from  Particles  and  Corresponding 
Power  Spectrum  Showing  Lorentzian  Line  Shape.  (The  Lorentzian 
halfwidth  implies  a  mean  diameter  of  3.4  Mm). 
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particle  size  distribution  mode  in  the  flow  is  not  much  greater  than  in  the 
original  distribution  (Fig.  4).  Due  to  time  constraints,  analytical  predic¬ 
tions  were  based  on  size  distributions  represented  by  three  size  groupings: 
0-1  um,  1-5  pm,  and  5-9  pm.  fable  1  is  n  typical  example  of  representative 
particle  sizes  for  each  grouping  and  their  relative  number  frequencies  for  a 
distribution  with  a  geometric  mean  of  2.0  pm  and  a  mode  of  0.5  pm.  This 
distribution  might  be  termed  "optimistic"  since  the  mean  is  smaller  than  indi¬ 
cated  by  the  size  diagnostic  measurements.  However,  as  will  be  discussed 
later,  the  resulting  calculated  spectra  are  relatively  insensitive  to  this 
sort  of  variation  in  the  size  distribution. 


TABLE  L 

PAR  T I  CL  F  SIZF  DISTRIBUTION 


Size  (pm) 

0.5 

3.5 

6.5 

Relative  Number 

BOO 

100 

5 

2.4  DATA  HANDLING  AND  REDUCTION 

Infrared  emission  and  laser  scattering  data  are  recorded  on  a  strip 
chart,  dual  pen  recorder  along  with  wavelength  event  marks.  This  information 
is  then  transcribed  to  digital  form  and  stored  on  a  computer  disc. 
Manipulation  and  correction  of  the  data  are  performed  by  computer  and  the 
corrected  output  is  then  graphed  by  a  computer-operated  plotter.  The  data  are 
primarily  corrected  for  the  optical  system  with  the  same  geometric  arrangement 
as  for  the  particle  flow.  The  measured  spectrum  of  the  calibration  source  is 
then  compared  to  the  "known"  spectrum  of  the  source  generated  by  computer  from 
the  Planck  equation:® 

UU,T)  =  ^5  exp(hc/XkT)-l  *  (16) 


8.  Reif,  F.,  Fundamentals  of  Statistical  and  Thermal  Physics  (McGraw-Hill, 
New  York,  1965) ,  pp.  3^8-388. 
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Differences  between  the  two  are  due  to  the  response  function  of  optical  system 
and  detector.  Once  the  response  function  was  determined,  it  was  stored  in  the 
computer  for  reduction  of  all  subsequent  data  taken  of  tieated  particle  flow. 
Various  tests  on  the  entire  reduction  system  indicate  that  the  amount  of  error 
introduced  by  the  digitization  and  computer  reduction  was  held  to  within 
1.0  percent  and  typically  runs  about  0.2  percent.  Since  the  signal-to-noise 
ratio  in  the  particle  data  runs  between  five  and  twenty  percent,  this  error 
introduced  by  processing  is  negligible.  As  mentioned,  the  spectrum  from  1.5 
to  5.5  urn  wavelengths  had  to  be  taken  in  three  segments  with  cutoff  filters 
primarily  to  avoid  difficulties  with  second  order  diffraction.  For  the  long 
wavelength  region  (3.6  to  5.5um)  a  different  grating  waa  also  required. 
Figures  8,  9,  and  10  illustrate  the  relative  instrument  response  curves  for 
the  three  wavelength  regions  respectively.  In  each  case  the  blackbody 
calibration  source  is  recorded  for  a  number  of  different  temperatures.  These 
spectra  are  then  divided  by  the  appropriate  blackbody  (Planck)  curve  to  pro¬ 
vide  the  instrument  response  curves  illustrated.  The  response  curves  from  the 
different  temperatures  are  averaged  together  to  improve  signal-to-noise  ratios 
before  being  used  to  reduce  data  from  the  particle  flow.  Since  the  radiation 
from  the  calibration  source  is  known  in  absolute  intensity  (not  just  relative 
spectral  intensity),  it  was  possible  to  obtain  absolute  instrument  response. 
(Note  that  ttv»  vertical  scales  in  Figures  8  through  10  are  in 

millivolts/milliwatts. )  In  the  calibration  procedure  it  was  not  physically 
possible  to  locate  the  calibration  source  at  the  focal  volume  of  the  monochro¬ 
mator  and  collecting  lens.  Instead,  an  aperture  was  placed  in  the  focal 
volume  and  the  calibration  source  was  located  80  mm  behind  the  aperture. 
Since  the  monochromator  was  operated  with  0.5  mm  x  20  mm  slits  and  2:1 
imaging,  the  aperture  used  for  calibration  was  0.25  mm  x  10  mm.  Therefore, 
the  measured  spectra  in  millivolts  was  divided  by  the  Planck  equation  (Watts) 
and  multiplied  by  the  solid  angle  of  the  radiation  times  the  area  of  the 

aperture.  It  should  be  noted  that  the  angular  dependence  has  been  removed 
from  the  common  form  of  Planck's  radiation  equation  (Eq.  16)  by  integrating 
the  solid  angle  over  a  hemisphere.  The  actual  angular  dependence  is 

cos  Q  where  9  is  tfye  angle  between  the  observer  and  the  normal  to  the 

radiating  surface  (calibration  source).  Hence,  the  integration  results  in  a 
factor  of  tt  in  Eq.  (16).  However,  due  to  the  small  aperture  used,  the 
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Figure  10. 


cos  8  term  can  be  accurately  approximated  by  1,  and  therefore  it  is  legitimate 
to  simply  multiply  by  the  subtended  solid  angle  rather  than  make  some  attempt 
at  integration. 
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SECTION  3 

RESULTS  AND  ANALYSIS 
3.1  PREDICTIVE  CALCULATIONS 

In  order  to  evaluate  the  data  obtained  from  thermal  emission,  it  was 
necessary  to  compare  the  data  with  theory.  The  most  convenient  way  to  do  this 
is  to  generate  calculated  spectra  using  the  measured  particle  size,  measured 
temperature,  and  established  values  for  the  index  of  refraction.  The  calcu¬ 
lated  spectra  were  then  compared  to  the  reduced  measurements.  Since  part  of 
the  initial  impetus  for  this  experiment  was  to  confirm  or  refute  a  novel  pre¬ 
diction  for  particle  radiation  (discussed  in  the  following  paragraphs),  the 
theoretical  aspects  are  presented  first. 

3.1.1  Modified  Planck  Eguationa 

from  basic  physics,  the  discussion  of  thermal  radiation  usually 
begins  with  an  enclosed  cavity  whose  walls  are  maintained  at  a  constant  tem¬ 
perature  so  that  with  a  sufficient  amount  of  time  the  radiation  field  inside 
the  cavity  comes  into  equilibrium  with  the  walls.  This  is  a  well-defined  and 
straightforward  problem  in  equilibrium  thermodynamics,  and  the  energy  density 
of  the  radiation  inside  the  cavity  is  given  by  the  Planck  equation9 


>>/  t\j  _  hw^du  w 

U(w,T)dw  =  — -  - 

2x  c  exp(hu/2xkT)  -  1 


(17) 


where  the  energy  per  unit  volume  and  unit  frequency,  U,  is  a  function  of  the 
frequency  (w)  and  teinperature  (T),  h  is  Planck's  constant,  k  is  Boltzmann's 
constant,  and  c  is  the  velocity  of  light.  It  is  fairly  common  to  generalize 
this  problem  by  considering  a  small  hole  or  aperture  in  the  cavity  through 
which  radiation  can  pass,  and  to  describe  the  absorption  and  emission  of  an 
"ideal"  object  (blpckbody)  as  being  represented  by  the  radiation 
absorption/emission  of  the  aperture.  The  radiation  flux  emitted  by  the  cavity 


9.  Ibid. 
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through  the  aperture  is  obtained  fro*  Eq.  (17)  by  including  a  factor  of 
c/4, and  can  be  rewritten  in  terns  of  the  wavelength,  as  Eq.  (16). 

In  Eq.  (17)  ihe  right  side  has  been  written  as  the  product  of  two  parts. 
The  first  part  describes  the  number  of  photons  per  mode  (eigenstate  in  the 
cavity)  and  the  second  part  comes  from  the  density  or  number  of  nodes 
available  in  the  cavity.  A  recent  paper*!  focused  on  the  degree  of  arbitrari¬ 
ness  with  which  this  density  of  modes  is  derived  and  suggested  that  it  be 
altered  to  include  the  situation  in  which  radiation  wavelengths  are  comparable 
to  cavity  dimensions.  The  modified  form  of  the  Planck  equation  was  thus 
derived  and  has  the  following  form 


MjU.T.d)  =  U(X,T)  {1  -  ( X2/2ir2d2) }  (18) 

where  is  the  new  radiance  function,  U  is  given  in  Eq.  (16),  and  d  is  the 
cavity  diameter.  Since  data  have  previously  been  scarce,  it  was  speculated 
that  this  formulation  could  be  applied  to  particle  thermal  radiation.  For 
diameters  less  than  one  am,  the  effect  of  thiB  new  term  is  very  significant. 
The  radiation  actually  terminates  instead  of  tailing  off  to  wavelengths  longer 
than  Sym. 


3.1.2  Hie  Theory  Calculations 

The  conventional  method  of  calculating  thermal  radiation  is  to 
multiply  the  standard  Planck  equation  by  an  emissivity  function.  In  general, 
the  emissivity  is  characteristic  of  the  material  in  question  and  is  a  way  of 
accounting  for  the  fact  that  a  given  material  may  absorb  or  emit  more  or  less 
strongly  in  one  wavelength  region  than  in  another,  in  comparison  to  the 
"ideal"  blackbody  described  by  the  Planck  equation.  The  emissivity  of  a 
material  is  related  to  the  index  of  refraction  (dielectric  constant)  of  the 
material  and  hence  is  measurable.  However,  in  the  situation  where  the  dimen¬ 
sions  of  the  material  object  are  comparable  to  the  wavelengths  of  interest  (as 


10.  Kittell,  C.,  Thermal  Physics  (Wiley,  New  York,  1969),  p.  256. 

11.  McGregor,  W.  K.,  up.  cit. 
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ia  the  case  with  thermal  radiation  from  pm-sized  particles),  geometric  optica 
is  no  longer  applicable,  and  some  attempt  must  be  made  to  solve  the  electro¬ 
magnetic  boundary  value  problems. 

The  solution  to  Maxwell's  equations  for  a  plane  electromagnetic  wave 

impinging  on  a  sphere  with  a  given  dielectric  constant  is  generally  referred 

to  as  Hie  Theory,  the  details  of  which  are  discussed  elsewhere. 12  Mie  Theory 

provides  a  method  to  calculate  the  total  extinction  cross  section  and  the 

total  scattering  cf oss  section  of  a  sphere  as  a  function  of  radius, 
wavelength,  and  index  of  refraction.  The  difference  between  the  amounts  of 
extinction  and  scattering  is  by  definition  the  amount  of  absorption. 
Therefore,  an  absorption  efficiency  factor  (Qabg  =  absorption  cross  section/ 
geometric  cross  section)  can  be  calculated.  Thermal  radiation  of  a  particle 
is  then  obtained  by  multiplying  the  efficiency  factor  (equal  to  the  emissivity 
via  Kirchoff's  Law)  by  the  Planck  blackbody  radiation  equation,  Eq.  (16): 


M2(A,T,d,n)  =  Qabs(M,n)  U(A,T) 


(19) 


where  n  is  the  bulk  index  of  refraction.  Calculations  of  Q.  used 

abs 

1.7  -  i  0.80  for  the  index  of  refraction  for  carbon  particles. For  larger 
particles  (1.0pm  diameter  and  above),  Qgba  is  fairly  flat  across  this  wave¬ 
length  range  so  that  the  resulting  particle  radiation  will  not  be  signifi¬ 
cantly  different  from  the  normal  Planck  radiation  spectrum,  Eq.  (16). 
However,  for  smaller  particles  this  is  not  the  case;  the  Planck  radiation 
curve  will  be  heavily  weighted  at  shorter  wavelengths. 


3.1.3  Theoretical  Predictions 

Computer  programs  have  incorporated  equations  17  and  18  to  provide 
graphs  of  the  calculated  thermal  radiation  spectra.  To  illustrate  the 


12.  Kerker,  M.,  Th^  Scattering  of  Light  and  Other  Electromagnetic  Radiation 

(Academic  Press,  New  York,  1969).  ’ 

13.  Pluchino,  A.  B,,  Goldberg,  5.  S,,  Dowling,  3.  M.,  and  Randall,  C.  M., 
3.  Appl.  Opt.  19,  3J70  (1980). 
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differences  involved,  Figure  11  compares  three  calculated  spectra  at  a  temp¬ 
erature  of  800  K.  The  solid  curve  corresponds  to  the  Planck  blackbody 
radiation  law  (Fq.  16).  The  lower  daahed  curve  corresponds  to  the  modified 
form  of  the  Planck  equation  (Eq.  18)  summed  over  the  particle  size  distribu¬ 
tion  shown  in  Figure  A.  To  show  this  spectrum  on  the  same  vertical  scale,  the 
values  had  to  be  nultiplied  by  a  factor  of  10.  The  upper  dotted  curve  is 
calculated  using  Mie  theory  (Eq.  19)  summed  over  the  same  particle  size 
distribution  and  using  an  index  of  refraction  of  1.7  -  i  0.8.  The  spectra 
reflect  the  interesting  features  previously  discussed  in  which  the  modified 
Planck  equation  predicts  a  definite  long  wavelength  cutoff  to  the  radiation, 
and  in  which  the  absorption  coefficient  from  Mie  theory  emphasizes  the  shorter 
wavelengths  compared  to  the  longer  wavelengths.  Both  of  these  effects  are 
dependent  on  small  particle  size.  It  should  be  expected  that,  aa  this  para¬ 
meter  is  increased,  both  broken  curves  will  eventually  converge  to  the  solid 
curve  since,  physically,  a  transition  is  being  made  from  a  small  particle  to  a 
macroscopic  object  for  which  the  Planck  equation  (Eq.  16)  is  known  to  be 
accurate.  Figure  12  illustrates  this  assertion  with  graphs  of  Eq's  (16), 
(18),  and  (19)  again,  except  that  the  particle  distribution  used  for  Eq's  (18) 
and  (19)  contains  the  larger  particle  sizes  shown  in  Table  I  at  a  temperature 
of  800  K. 

To  compute  the  thermal  radiation  from  a  size  distribution,  each  par¬ 
ticle  size  is  taken  separately  and  the  absolute  radiation  spectrum  is  calcu¬ 
lated  for  that  size,  Thi3  spectrum  is  then  multiplied  by  the  surface  area 
and  the  relative  number  frequency  of  the  particle  in  the  distribution.  These 
separate  spectra  are  then  added  together  for  each  wavelength  value  and  the 
total  is  divided  by  the  total  surface  area  of  the  distribution. 
Algebraically,  Figures  11  and  12  display  the  quantity: 
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where  is  the  appropriate  radiance  function  (Eq's  18  or  19),  the  index,  i, 
refers  to  the  particle  size,  S  is  the  surface  area  of  the  particle,  and  N  is 


* 


I  I  I 
* 
o 


— 1 

toil. 

O^Q 
UJ  ^  IjJ 

"38 

2  0-  5 


*\ 


\  *•-. 

\  *•. 
\ 

\ 


(U0!P0J3|$  .  _UiO  .  UlTf/  SUDM  )  A1ISN31NI 


Calculated  Spectra  of  Infrared  Emission  from:  a  Blackbody  (Solid  Curve),  Particles 
Using  Hie  Theory  (Dotted  Curve),  and  Particles  Using  the  Modified  Planck  Equation 
(Dashed  Curve  Multiplied  by  10  to  Graph  on  the  Same  Scale).  Particle  size  distribution 
is  shown  in  Figure  4,  computed  at  800  K. 


Calculated  Spectra  of  Infrared  Emission  from:  a  Blackbody  (Solid  Curve),  Particles 
Using  Hie  Theory  (Dotted  Curve),  and  Particles  Using  the  Modified  Planck  Equation 
(Dashed  Curve).  Particle  size  distribution  is  shown  in  Table  1,  computed  at  800  K. 


its  number  frequency  in  the  size  distribution.  The  reason  the  modified  Planck 
formulation  predicted  a  much  lower  intensity  in  figure  11  is  now  easily  seen, 
for  the  very  email  particle  sizes  shown  in  figure  4,  the  factor 
in  Eq.  (18)  diminishes  or  completely  cuts  off  any  contribution  to  the  total 
radiation  for  more  than  90  percent  (by  number)  of  the  particles  at  wavelengths 
of  2.5  pm  and  larger.  Therefore,  compared  to  Eq.  (19),  the  total  radiation 
intensity  will  be  substantially  less  for  the  same  particle  sample.  As  the 
particle  size  becomes  larger,  however,  this  factor  becomes  less  important,  and 
as  seen  in  figure  12,  the  modified  Planck  formulation  predicts  a  similar  peak 
intensity  for  the  large  particle  size  distribution  given  in  Table  1. 

3.2  EXPERIMENTAL  RESULTS 

Several  sets  of  spectra  were  taken  for  each  wavelength  region.  In  each 
set  the  spectra  were  taken  consecutively  in  time  so  that  the  experimental 
parameters  of  the  run,  flow  velocity,  temperature,  and  particle  charac¬ 
teristics  were  held  as  constant  as  possible.  The  reduced  data  shown  in 
figure  13  are  an  example  of  a  aeries  of  spectra  taken  consecutively  showing 
the  "shot -to-shot'1  reproducibility.  These  sets  of  spectra  are  then 
numerically  averaged  to  reduce  random  noise.  Averaged  sets  from  the  different 
wavelength  regions  are  then  matched  on  the  basis  of  the  measured  flow  tem¬ 
perature  and  slope  at  the  point  of  juncture,  figure  14  is  an  example  of  com¬ 
posite  data  in  which  averaged  spectra  from  different  wavelength  regions  have 
been  graphed  on  the  same  wavelength  scale.  for  convenience  in  comparison, 
predictive  calculations  shown  previously  are  included.  The  dotted  curve  is 
the  Mie  theory  prediction  while  the  dashed  curve  is  the  modified  Planck 
equation  curve.  The  two  calculated  curves  have  been  normalized  to  the  same 
height  to  match  the  data,  and  both  calculations  employ  the  same  particle  size 
distribution  shown  in  Table  I  for  a  temperature  of  820  K  at  which  the  data 
were  taken.  As  predicted  in  figure  12,  the  two  calculated  spectra  are  very 
similar.  Both  are  in  good  agreement  with  the  measured  spectrum. 

for  a  large  particle  size  distribution,  a  distinction  between  the  Mie 
theory  approach  and  the  modified  Planck  equation  is  difficult  to  achieve  on 


consecutively  during  a  single 
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the  uasis  of  relative  spectral  shape,  given  the  noise  level  of  the  data  and 
the  degree  of  error  in  the  experimental  parameters.  Another  possibility  was 
differentiating  the  predictive  formulations  on  the  basis  of  absolute  intensity 
of  radiation.  Since  the  optical  system  was  calibrated  against  a  standard 
blackbody  source,  it  was  possible  to  determine  the  absolute  power  levels  of 
the  infrared  data  from  the  particles  as  described  below.  Significant 
variation  from  different  data  sets  is  accounted  for  by  differences  in  flow 
temperature  and  particle  density.  Nearly  all  the  data  sets  yield  peak  inten¬ 
sity  levels  between  0.1  and  1.0  mM/iim.  Since  laser  attenuation  meaauremsnts 
were  made  during  some  of  the  experimental  runs  the  quantity: 


-ln(I/I0)/L  =  £  Nj.j' 


is  easily  determined,  where  1  and  I  are  the  initial  and  attenuated  laser  beam 
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intensities,  L  is  the  pathlength  through  the  flow,  and  and  o  ^  are  the 
number  density  and  extinction  cross  section  for  particles  of  size  i.  From 
Eq.  (21)  it  is  straightforward  to  estimate  the  total  geometric  surface 
area  (S)  of  all  the  particles  present  in  the  focal  volume  (V)  of  the  scanning 
monochromator: 


s  =  *E*i'V!Evi 


where  a  is  the  geometric  cross  section  of  the  i  n  size  particle 

«  *2  in 

(o  j  =  xri  ),  and  from  Hie  theory,  o  i  “20^  Using  this  surface  area,  the 
measured  focal  volume,  and  the  subtended  solid  angle  of  the  collection  optica, 
the  Mie  theory  approach  predicts  radiated  intensities  of  around  0.1  mN/unit 
wavelength,  for  both  the  large  and  small  size  distributions.  The  modified 
Planck  equation  predicts  values  20  times  lower  (see  Figure  11)  for  the  small 
size  distribution,  bqt  values  only  about  25  percent  lower  for  the  large  size 
distribution  (see  Figpre  12).  Consequently,  both  the  Hie  theory  and  modified 
Planck  equation  formulations  predict  radiation  inteneities  for  a  large  size 
distribution,  which  are  within  experimental  error  of  the  measured  values. 


M 


36 


This  approach  only  provides  a  means  of  discrimination  in  the  small  particle 
size  regime  where  differences  in  spectral  shape  should  also  be  significant. 
Therefore,  the  immediate  result  is  that  the  experimental  data  do  not  permit  a 
clear  distinction  between  the  two  theoretical  formulations  presented  here  due 
to  the  large  (>  1.0  urn)  average  particle  size  in  the  two-phase  flow. 

3.3  ANALYSIS 

Three  basic  parameters  or  experimental  variables  are  essential  to  predic¬ 
tive  calculations  of  thermal  radiation  from  particulate  material.  They  are 
size,  emissivity  (or  index  of  refraction),  and  temperature.  For  the  experi¬ 
ments  described  here,  each  of  these  variables  was  measured  independently  to 
generate  the  calculated  spectra  which  generally  agree  with  the  observed 
radiation  spectra.  The  question  now  considered  is  this.  Can  any  information 
about  these  variables  be  derived  from  the  observed  spectra  for  situations 
where  prior  knowledge  is  not  easily  obtained?  More  simply,  what  is  the 
diagnostic  value  of  the  thermal  emission?  This  question  can  be  explored  by 
using  the  theoretical  calculations  which  were  seen  to  be  consistent  with  the 
data.  Figure  13  shows  the  dependence  of  the  thermal  emission  from  0.5,  1.0, 
2.0,  and  4.0  pm  diameter  carbon  particles  at  a  temperature  of  800  K  using  the 
formulation  based  on  Mis  theory.  These  spectra  for  different  single  size  par¬ 
ticles  have  all  been  normalized  to  approximately  the  same  height  to  show  the 
relative  spectral  distributions.  Differences  in  absolute  radiation  inten¬ 
sities  were  less  than  a  factor  of  two.  From  Figure  15,  it  is  suggested  that 
particles  larger  than  about  a  micrometer  in  diameter  have  very  similar 
spectra.  Indeed,  for  applications  where  the  emitted  wavelength  is  larger  than 
or  equivalent  to  the  particle  diameter,  this  result  is  even  independent  of  the 
theoretical  formulatipn,  since  the  transition  to  a  macroscopic  object  ia  being 
made  and  the  radiated  spectrum  is  dominated  by  the  Planck  equation. 
Therefore,  the  size  parameter  appears  to  be  relatively  unimportant  for  ^pli¬ 
cations  involving  larger  sizes  of  particles  at  high  temperatures. 

The  next  parameter  examined  is  the  emissivity,  or  in  the  Mie  theory 
formulation,  the  index  of  refraction  of  the  particulate  material.  The  impor¬ 
tant  aspect  here  is  the  interdependence  of  the  optical  properties  of  materials 
with  the  particle  size.  While  this  is  a  fundamental  question,  subject  to  much 


interest  and  active  investigation,  the  data  collected  for  carbon  particles 
suggests  that  no  dramatic  changes  occur  in  the  index  of  refraction  for  par¬ 
ticle  sizes  on  the  order  of  2  pm.  Figure  16  displays  calculated  spectra 
(using  Mie  theory)  for  a  single  particle  size  of  2.0  um,  using  values  of  1.0, 
1.7,  and  3.0  for  the  real  part  of  the  index  of  refraction.  (The  imaginary 
part  was  -0.8.)  The  three  spectra  are  graphed  on  an  absolute  scale  and  indi¬ 
cate  that  again  relatively  little  effect  is  made  by  changing  the  index  of 
refraction  (real  part)  by  factors  of  almost  two.  Similar  changes  in  the  imag¬ 
inary  part  of  the  index  of  refraction  had  a  much  greater  effect  on  the  abso¬ 
lute  intensity  of  radiation,  but,  when  normalized,  had  negligible  effect  on 
the  relative  spectral  distribution.  Barring  materials  that  have  dramatically 
different  optical  constants  (such  as  near  zero  imaginary  parts  of  the  index  of 
refraction),  and  considering  particles  in  the  aize  range  of  a  few  micrometers 
or  more,  one  could  generalize  these  results  to  aay  that  the  thermal  radiation 
spectral  distribution  is  not  a  strong  function  of  material  composition. 

The  last  parameter  is  the  temperature  of  the  two-phase  flow.  Figure  17 
shows  another  composite  graph  of  reduced  data  overlaid  with  three  theoretical 
curves.  The  three  spectra  were  calculated  using  the  same  particle  aize 
distribution  (Table  I),  but  at  temperatures  of  800,  850,  and  900  K.  The  data 
themselves  were  recorded  with  the  measured  temperature  in  the  two-phase  flow 
of  840  ±  18  K.  It  can  be  seen  that  the  theoretical  curve  of  850  K  appears  to 
be  a  best  fit  to  thq  data,  while  the  other  two  curves  car.  be  qualitatively 
rejected  as  marginal  fits.  ThiB  is  an  approximate  way  of  establishing  the 
accuracy  with  which  the  temperature  can  be  derived  from  the  data.  In  the 
absence  of  significant  dependence  on  the  other  two  parameters,  size  and  index 
of  refraction,  the  temperature  can  be  "deduced"  from  the  data  within  about 
±  6  percent  (50/850),  This  result  invites  speculation  as  to  whether  such 
spectra  could  provide  temperature  measurements  in  applications  such  as  rocket 
plumes.  Figure  18  phows  calculated  spectra  for  a  2.0  pm  carbon  particle 
exactly  as  in  Figure  17,  but  with  the  temperatures  elevated  to  around  3,100  K. 
The  graphs  in  Figure  18  suggest  that  with  similar  signal-to-noise  ratios,  the 
thermal  radiation  frpm  a  two-phase  flow  could  provide  measurement  of  these 
elevated  tenperatures  with  better  than  10  percent  accuracy. 


Figure  17.  Co«>osite  Data  Graphed  on  a  Single  Wavelength  Axis  Compared  to  Three  Calculated 

Spectra  U8ing  Hie  Theory.  The  teaperature  of  the  data  ia  840  K,  irfiile  the  calculations 
were  wde  for  800,  850,  and  900  K.  Calculations  used  the  particle  size  distribution 
shown  in  Table  1. 


SECTION  4 
CONCLUSIONS 


In  terms  of  tha  ultimate  applications  of  this  work,  it  appears  that  the 
two  fundamental  theoretical  approaches  for  determining  particle  emiasivitiea 
are  not  a  crucial  factor  in  predicting  infrared  emission  as  long  as  the  con¬ 
ditions  of  large  particles  size  (>  1  urn),  and  high  temperatures  (>  800  K)  are 
met.  In  this  regime  theoretical  calculations  of  the  thermal  radiation  will  be 
dominated  by  the  Planck  equation  for  a  large  range  of  particle  sizes  and  a 
large  range  of  matf rials.  This  further  implies  that  the  measurement  of 
infrared  emission  may  be  used  as  a  diagnostic  probe  for  the  temperature  of  the 
flow.  For  particular  application  to  rocket  exhausts,  calculations  indicated 
that  at  higher  temperatures  (3000  K)  an  accuracy  of  t  10S  or  better  may  be 
achieved.  These  results  were  obtained  from,  and  therefore  are  directly  appli¬ 
cable  only  to  two-phase  flows  which  are  optically  thin.  Complications  due  to 
radiative  transfer  have  not  been  considered.  In  ultimate  application  to  opti¬ 
cally  thick  flows  this  last  consideration  together  with  large  spatial  gra¬ 
dients  may  be  the  limiting  factors. 
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